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ABSTRACT 

We present observations with the Chandra X-ray Observatory of the pulsar wind nebula (PWN) powered by 
the energetic young pulsar B1509-58. These data confirm the complicated morphology of the system indicated 
by previous observations, and in addition reveal several new components to the nebula. The overall PWN shows 
a clear symmetry axis oriented at a position angle 150° ± 5° (north through east), which we argue corresponds to 
the pulsar spin axis. We show that a previously identified radio feature matches well with the overall extent of the 
X-ray PWN, and propose the former as the long-sought radio nebula powered by the pulsar We further identify a 
bright collimated feature, at least 4' long, lying along the nebula's main symmetry axis; we interpret this feature 
as a physical outflow from the pulsar, and infer a velocity for this jet > 0.2c. The lack of any observed counter-jet 
implies that the pulsar spin axis is inclined at ^ 30° to the line-of-sight, contrary to previous estimates made from 
lower-resolution data. We also identify a variety of compact features close to the pulsar. A pair of semi-circular 
X-ray arcs lie 17" and 30" to the north of the pulsar; the latter arc shows a highly-polarized radio counterpart. 
We show that these features can be interpreted as ion-compression wisps in a particle-dominated equatorial flow, 
and use their properties to infer a ratio of electromagnetic to particle energy in pairs at the wind shock a ^ 0.005, 
similar to that seen in the Crab Nebula. We further identify several compact knots seen very close to the pulsar; 
we use these to infer a < 0.003 at a separation from the pulsar of 0. 1 pc. 

Subject headings: ISM: individual: (G320.4-1.2) — ISM: jets and outflows — pulsars: individual (B1509-58) 
— stars: neutron — supernova remnants — X-rays: ISM 
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1. INTRODUCTION 

Pulsars and supernova remnants (SNRs) are both believed to 
be formed in supernova explosions. However, there are still 
fewer than 20 cases in which convincing associations between 
a pulsar and a SNR have been established. These few systems 
provide important information on the properties of pulsars and 
SNRs, and on the relationship and interaction between them. In 
particular, in cases where a pulsar is still inside its associated 
SNR, the high-pressure environment can confine and shock the 
pulsar's relativistic wind. This generates synchrotron emission, 
producing an observable pulsar wind nebula (PWN) which can 
be used to trace the energy flow away from the pulsar 

For almost 15 years, the Crab Nebula and the Vela SNR 
were the only SNRs known to be associated with pulsars. A 
third association only emerged when a 150-ms X-ray, radio 
and 7-ray pulsar, PSR B1509-5 8, was discovered within the 
SNR G320.4-1.2 (MSH 15-52) (Reward & Harnden Jr. 1982 



Manchester, Tuohy, & D' Amico 1982| ; |Ulmer et al. 1993D . The 
pulsar's spin-parameters can be used to infer a characteristic 
age Tc = 1700 yr, a spin-down luminosity £ = 1 .8 x lO-'^ erg s"' 
and a surface magnetic field B„ = 1.5 x 10'-^ G (K aspi et al. 
1994), making it one of the youngest, most energetic and 



highest-field pulsars known. 

SNR G320.4-1.2 has an unusual radio appearance, consist- 
ing of two well-separated radio regions (Fig, [ll; Gaensler et al. 
1999, hereafter G99). While the more southerly of these re- 
gions approximates a partial shell, the northern component is 
bright, centrally concentrated, and has an unusual ring of radio 
clumps near its center. While it had been previously suggested 
that G320.4-1.2 consisted of two or three separate SNRs, G99 
have used Hi absorption measurements to confirm that this 
source is indeed a single SNR at a distance of 5.2± 1.4kpc. As- 
suming standard parameters for the interstellar medium (ISM) 
and for the supernova explosion, an age of ~6-20 kyr is derived 
for the SNR, an order of magni tude larger than the 1700 yr im- 
plied by the pulsar spin-down ( Seward et al. 1983[ ). One pos- 
sible explanation for this discrepancy is that the p ulsar is much 
olde r than it seems (|Blandford & Romani 1988|; Gvaramadze 
2001), but this requires an unusual spin-do wn history of whic h 
there is no evidence from long-term timing ( Kaspi et al. 1994| ). 



A more likely explanation is that the SNR has expanded rapidly 
into a low-density cavity ( Seward etal. 198^ ). This model can 



also explain the unusual SNR morphology, the offset of the pul- 
sar from the SNR's appare nt center, and the f aintness of any 
PWN at radio wavelengths (|Bhattacharya 1990|; G99). This ar- 
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gument is supported by recent observations of H I emission in 
the region (Dubner et al. 2002). 

X-ray emission from this system, also displayed in Figure |l], 
has revealed a different but similarly complicated picture. Apart 
from a pulsed point-source c orresponding to the pulsar itself 
(Seward & Harnden Jr. 1982), X-ray emission from the system 
is dominated by an elongated non-thermal source centered on 
the pulsar, and an extended thermal source coincident with the 
northern radi o component of the SNR and with the optic al neb- 



ula RCW 89 dSeward et al. 1983| ; [Trussoni et al. 1996| ). Faint 



diffuse X-rays extend over the entire extent of the radio SNR, 
and potenti ally correspond to th ermal emission from the SNR 
blast wave (Trussoni et al. 1996). 



The central non-thermal so urce has been interp reted as a 
PWN powe red by the pulsar (^ward et al. 1984|; T russoni 
et al. 1996). As mentioned above, no radio emission from 
this PWN has been seen, presumably bec ause of the low den- 
sity into which the pulsar wind expands ( [Bhattacharya 1990 ; 
G99). High-resolution observations with ROSAT show the X- 
ray PW N to be significantly collim ated to the SE of the pulsar 
(Fig. |l|; [Greiveldinger et al. 1995 ). Several authors have in- 
terprete d this region as corresponding to a jet prod uced by the 
pulsar ( [Tamura et al. 1996 ; Brazier & Becker 1997). 

Using ROSAT PSPC data, Greiveldinger et al. (1995) have 
analyzed the morphology of X-ray emission close to the pul- 
sar, and concluded that the PWN has two components: a large 
region extending over many arcminutes, and a compact disc of 
emission of radius 30" -45" centered on the pulsar Brazier 
& Becker (1997) considered data taken with the ROSAT HRl. 
They found no evidence for the two-component PWN proposed 
by Greiveldinger et al. (1995), but pointed out that the X- 
ray emission around the pulsar had a cross-shaped morphology. 
They interpreted this as indicating that the pulsar powers an 
equatorial torus, as is also seen in X-rays for the Crab pulsar 
(Aschenbach & Brinkmann 1975; [Hester et al. 1995 ). This 
morphology would imply that the jet (and presumably the pul- 
sar spin axis) lies largely in the plane of the sky, inclined at 
~ 70° from the line-of-sight. 

An interpretation for the thermal X-ray source ^ 6' to the 
north of the pulsar is unclear It has been repeatedly proposed 
that this region represents an interaction between a collimated 
outfl ow from the pulsar and the surrounding SNR (Sew ard et al. 
1983; [Manchester & Durdin 1983t [Manchester 1987t B razier 
& Becker 1997). Tamura et al. (1996) showed evidence for 
a bridge of non-thermal emission connecting the pulsar to the 
thermal X-ray emission to the north, but it remained uncertain 
as to whether this was the direct counterpart to the SE jet-like 
feature. At high spatial resolution, a collection of X-ray clumps 
is seen at the center of the northern X-ray region, the morpholo- 
gies and positions of which closely match thos e of clumps seen 
at radio wavelengths (Brazier & Becker 1997; G99). It has not 
been established whether t he proposed pulsar outflow powers 
the entire RCW 89 region (Tamura et al. 1996) or just the cen- 
tral collection of clumps (G99). But in either case, the evidence 
for interaction provides a convincing argument that the pulsar 
and SNR are physically associated. 

Clearly PSR B1509-58 and SNR G320.4-1.2 provide an im- 
portant example of the interaction between a pulsar and its 
environment. Accordingly, we have obtained sensitive high- 
resolution observations of this system with the Chandra X-ray 



Observatory. These data provide the opportunity to clarify the 
nature of the emission close to the pulsar, to determine the over- 
all structure of the PWN, and to further investigate the apparent 
interaction between the pulsar and its surrounding SNR. In §^ 
we describe our observations, in §|| we describe the resulting 
images and spectra, and in §Q we discuss and interpret these 
results. 

2. OBSERVATIONS AND ANALYSIS 

G320.4-1.2 was observed on 2000 Aug 14 (ObsID 754) with 
the Advanced CCD Imaging Spectrometer (ACIS) abo ard the 
Chandra X-ray Observatory (Weisskopf et al. 2000b). ACIS 
is an array of ten 1024x 1024-pixel CCDs fabricated by MIT 
Lincoln Laboratory. These X-ray sensitive devices have large 
detection efficiency (10-90%) and moderate energy resolution 
(10-50) over a 0.2-10.0 keV passband. Coupled with the 10-m 
optics of Chandra, the CCDs have a scale of 0."492 per pixel, 
well-matched to the on-axis point-spread function (FWHM ^ 
1"). ACIS and its calibration are discussed in detail by Burke 
etal. (1997) and Bautzetal. (1998). 

A single exposure of length 20 ks was made in the stan- 
dard "Timed Exposure" mode. The data discussed here were 
recorded from the four CCDs of the AClS-1 array. After stan- 
dard processing had been carried out at the Chandra Science 
Center (ASCDS version number R4CU5UPD14.1), we ana- 
lyzed the resultant events list using ClAO v2.1.3. We first 
applied the observation-specific list of bad pixels to the data, 
and then corrected for the ^1" offset present in ACIS -I obser- 
vations". After removing small intervals in which data were 
not recorded, the final exposure time for this observation was 
19039 s. 

The effective area of the detector is a function of both the 
energy of incident photons and their position on the sky. To 
generate an exposure-corrected image, we computed exposure 
maps at five energies: 0.5, 1.5, 3.0, 5.0 and 7.0 keV. We also 
computed images of the counts per pixel in the five separate en- 
ergy bands 0.3-0.8, 0.8-2.0, 2.0-4.0, 4.0-6.0 and 6.0-8.0 keV. 
For each energy band, we divided the map of raw counts by 
the exposure map corresponding to that band, and then normal- 
ized by the exposure time to produce a final image with units of 
photons cm"^ s"' . 

For compact sources of emission (< 5" in extent), spec- 
tra were extracted from the data using the psextract script 
within ClAO, using a small extraction region immediately sur- 
rounding the source of interest. For extended sources we made 
use of vl.08 of the calcrmf and calcarf routines'^ pro- 
vided by Alexey Vikhlinin, which generate the response ma- 
trix and effective area for a region of the CCD by computing 
an appropriately weighted mean over small sub-regions. In all 
cases, the resulting spectra were rebinned to ensure a minimum 
of 15 counts per spectral channel. The background spectrum 
for each region was measured by considering an annular region 
surrounding the source of interest. Subsequent spectral analysis 
was carried out using XSPEC v 1 1 .0. 1 . 

3. RESULTS 

3.1. Imaging 

An exposure-corrected and smoothed image in the energy 
range 0.3-8.0 keV, encompassingthe entire field-of-view of the 
ACIS -I array, is shown in Figure pi The central region is shown 



nttp 


//asc 


harvard 


edu/mta/ASPECT/f ix_of f set . cgi 




nttp 


//asc 


harvard 


edu/cont-soft/ software /calcrmf .1.08. html 



GAENSLER ET AL 



3 



on a logarithmic scale in the subpanel at lower-right. Count 
rates and surface brightnesses for various features (marked in 
Fig. H and discussed below) are given in Table |l]. 

Mai-ked in Figure || as feature A, PSR B1509-58 is clearly 
detected as a bright source near the center of the field-of-view. 
The extent of this source is consistent with the point-spread 
function of the telescope. The pulsar's X-ray position is mea- 
sured to be (J2000) RA 15^13'"55"64, Dec -59°08'09"2, with 
an uncertainty of ±0."5 in each coordinate. This differs by sev- 
eral ar csec from previous X - ray measurements of the pulsar po- 
sition (Seward et al. 1984; Brazier & Becker 1997), but is in 



excellent agre ement with positi ons measured from radio timing 
and imaging ( |Kaspi et al. 1994| ; G99). 

Because the pulsar is so bright, a significant number of 
counts from it strike the detector during the 41 ms per 3.2- 
sec frame during which the CCD is read out. This produces a 
faint trail of emission, marked as feature B, which runs through 
the pulsar and parallel to the north/south edge of the array (i.e. 
along the read-out direction). 

The morphology of the surrounding X-ray emission as im- 
aged with Chandra is consistent with that seen by ROSAT and 
ASCA at lower spatial resolution. In particular, we confirm 
t he diffuse nebula seen extending ^ 6' NW/SE of the pulsar 
( Greiveldinger et al. 1995|; iTrussoni et al. 1996|; T amura et al. 
1996), and the colle ction of X-ray clumps c oincident with the 
Ha nebula RCW 89 ( [Brazier & Becker 1997| ). Most of the neb- 
ula is elongated along a position angle 150° ±5° (N through 
E). In all future discussion, we consider this orientation to de- 
fine the main axis of the system. 

To the SE of the pulsar and superimposed on the diffuse neb- 
ular emission is a jet-like X-ray feature, denoted as feature C in 
Figure |[ Close to the pulsar, feature C is aligned with the main 
nebular axis, but at a distance ^ 2' from the pulsar it curves 
around to the east, before fading abruptly at a separation of ^ 4' 
from the pulsar. This feature is resolved and is approximately 
10" wide. 

In Figure H we show a comparison between our Chandra data 
and 1 .4-GHz radio observations taken with the Australia Tele- 
scope Compact Array (ATCA; G99). There is a good match 
between the morphologies of the X-ray PWN and the radio 
"tongue" reported in previous observations (M anchester & Dur- 
din 1983; G99), the radio extent being somewhat larger than 
the X-ray nebula. Furthermore, it is clear that feature C corre- 
sponds to a trough of reduced radio emission, confirming a sim- 
ilar correspondence seen when comparing ATCA and ROSAT 
observations (G99). 

We see no obvious counterpart in Figure ||to feature C on the 
other side of the pulsar. Rather, to the NW of the pulsar there is 
an elongated region of reduced emission, beginning 3' from the 
pulsar and extending for 3' along the nebular axis. We refer to 
this as feature D. The width of this fainter region is a function 
of distance from the pulsar: at its end nearest the pulsar, this 
feature is ~ 25" across, but it widens to 75" at its NW end. 

Approximately 30" to the north of the pulsar is feature E, a 
semi-circular arc of emission, approximately centered on the 
pulsar. This arc is resolved by Chandra, with approximate 
width 15". The arc is brightest immediately to the north of the 
pulsar, fading slowly to the east and more rapidly to the west. 
Along the main symmetry axis defined above, the pulsar is sep- 
arated from the inner edge of this arc by 30", and from the outer 
edge by 45". 

We have re-examined the 1.4- and 5-GHz images presented 
by G99 to look for a radio counterpait to feature E. In total 



intensity, this region is confused by emission from RCW 89 
(see Fig. ||), and we can identify no specific radio feature which 
might be associated. However, we have re-analyzed the 5-GHz 
polarization data of G99, deconvolving the data with a new 
maximum entropy algorithm sp ecifically designed for polari - 
metric observations (PMOSMEM; [Sault, Bock, & Duncan 1999| ). 
In Figure ^ we compare the resulting distribution of linear po- 
larization to the X-ray emission in the same region. While there 
is generally a high level of polarization from the area, there is 
a clear enhancement of the linearly polarized signal in an arc 
^ 30" to the north of the pulsar, whose extent and morphology 
match closely that of the X-ray emission from feature E. 

Returning to Figure ^ it can be seen that the perimeter of 
the overall diffuse nebula is reasonably well-defined. This is 
particularly the case to the west of the pulsar, where feature F, 
a sharp-edged protuberance, extends ~ 2' to the WNW of the 
main body of the nebula. A similar, but smaller and less well- 
defined feature may exist on the opposite side of the nebular 
symmetry axis. 

Figure ^ shows the X-ray emission immediately surrounding 
the pulsar at the full resolution of the data. In increasing order 
of distance from the pulsar (marked as A), we point out various 
features of interest: 

1 . Approximately 3" to the NW of the pulsar, a small clump 
of emission ^ 3" across. This feature is too far from 
the pulsar and too irregular in appearance to be due to 
asy mmetries in the w ings of the point-spread function 
(cf. Pavlov et al. 2001 ). The count rate from this source is 



approximately constant over the observation, so it cannot 
result from pulsar photons assigned to the wrong loca- 
tion on the sky (as might result from errors in the aspect 
solution). 

2. Lying 8" to the SE of the pulsar along the main nebular 
axis, another small clump of emission. 

3. Approximately 10" to the north of the pulsar, a circular 
clump of diameter 3". 

4. Immediately to the NW of feature 3, a similarly-sized 
clump. 

5. Superimposed on feature 4, a faint circular arc of emis- 
sion. This arc is approximately centered on the pulsar, 
with a radius of 17" and an angle subtended at the pulsar 
of ~ 110°. The midpoint of this feature lies along the 
main nebular axis. 

6. A faint unresolved source, positionally coincident with 
the 06.5III star Muzzio 10 (Prsatti & Muzzio 1980|; Se- 
ward et al. 1983). 

3.2. Spectroscopy 

The gross spectral characteristics of the system can be de- 
rived from Figure ^ where we show the field surrounding 
SNR G320.4-1.2 in four separate energy sub-bands. This set 
of images demonstrates that the northern ring of X-ray clumps 
coincident with RCW 89 have no detectable emission above 
^4 keV, while the PWN has significantly harder emission. The 
pulsar is clearly the hardest source in the field, although this 
may largely be due to pile-up in its spectrum (see further discus- 
sion below). These broad conclusions confirm the spectral de- 
compositions m ade from Einstein and ASCA observations ( Se- 
wardetal. 1983; [Tamura et al. 1996^ 
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At the higher spatial resolution of the Chandra data, some 
new spectral features become apparent. The 2.0-^.0 keV and 
4.0-6.0 keV images in Figure g| demonstrate that the PWN has 
an extended component coincident with the RCW 89 clumps. 
Furthermore, there is the suggestion from Figure ^ that the jet 
and outer arc (features C and E respectively) have harder spec- 
tra than the overall PWN. 

Spectral fits show that the point-source seen at 
RA 15'>13'"4P, Dec -59°11'45" is heavily absorbed, with 
Nh ^ 3 X 10^^ cm"^ for simple thermal and non-thermal mod- 
els. It therefore most likely represents an unrelated background 
source. 

3.2.1. PSR Bl 509-58 

The pulsar itself has a high X-ray flux and is expected to 
suffer from significant pile-up, in which multiple events arrive 
on a given CCD pixel during a single frame. Indeed the spec- 
trum of the pulsar has an excess of hard photons (presumably 
due to several lower-energy photons being recorded as a single 
event), and cannot be fit by any simple power-law model. As- 
suming a foreground absorbing column Nh ~ 1 x 10^^ cm"^, 
a power-law spectrum with photon index F ~ 1.4, and an un- 
absorbed flux ^ 6 x 10"'^ erg cm~^ s"' in the energy range 



0.1-2.4 keV ( 3reiveldinger et al. 1995 ), the Chandra proposal 
planning software'^ predicts >50% pile-up, with a resulting de- 
tected count rate 0.16 cts s"' over the energy range 0.5-10 keV. 
This in reasonable agreement with the value measured for the 
pulsar in Table |l]. 

3.2.2. The Pulsar Wind Nebula 

We can determine the spectrum for the diffuse nebula sur- 
rounding the pulsar by extracting photon energies from the an- 
nular region shown in Figure ^ Measuring the spectrum for 
an extended source such as this is difficult because of radiation 
damage to the front-illuminated CCDs,''*, which has caused the 
response of each CCD to be a function of distance from the 
read-out nodes. To try and mitigate this effect, we extracted 
separate spectra for each of the four CCDs on which this dif- 
fuse emission falls, and generated response matrices and effec- 
tive area curves for each CCD separately. We then fit these 
four spectra simultaneously over the entire usable energy range 
(0.5-10.0 keV), using a common value for all fit parameters ex- 
cept the normalization of each spectrum. These spectra, which 
between them represent a total of ^^50 000 photons, are shown 
in Figure ^. While no significant spectral features are seen, 
some systematic residuals are still present, which we attribute 
to gain mismatches resulting from the aforementioned radia- 
tion damage. The count rate from the diffuse nebula is very 
high (^3 cts s"'), but this emission is spread over a very large 
area so that pile-up is negligible. 

We find that these data are well fitted by an absorbed power 
law. As listed in Table ^, the best-fit spectral parameters are an 
absorbing column Nh = (9.5 ± 0.3) x 10^' cm"^ and a photon 
index r = 2.05 ±0.04. 

We next consider the spectra of each of regions C and E and 
of features 1-5. All these sources have a surface brightness far 
too low to produce pile-up in the detectors. We fit absorbed 
power laws to the data in each region over 0.5-10.0 keV, with 
the resulting spectral parameters listed in Table ^. With the ex- 
ception of feature 2 (which has very few counts and correspond- 
ingly uncertain spectral parameters), the fitted values of Nh are 



all consistent with each other and with that for the diffuse PWN, 
while the resulting photon indices range between F « 1.3 and 
F w 1.7. These photon indices are all significantly flatter than 
that determined for the diffuse PWN. To more tightly constrain 
the photon index in these regions, we assume that they all have 
the same absorbing column as that determined for the diffuse 
nebula above, and thus fix Nh to 9.5 x 10^' cm~^ correspond- 
ing to the power law fit to that source. The resulting spectral 
fits confirm the flatter spectra for these regions when compared 
to the overall nebula. 

Finally, we crudely fit a Raymond-Smith spectrum to the 
emission from feature 6, coincident with the star Muzzio 10. 
The approximate spectral parameters are listed in Table ^ 

3.2.3. RCW 89 

The spectrum for the RCW 89 region is considerably more 
complex than for the PWN. In Figured we show a spectrum for 
the region shown in Figure ^, enclosing the brightest clumps 
in RCW 89 (sourc es NO, Nl, N2, N3, N5, N6, N7 and N8 in 
the designation of Brazier & Becker 1997). This spectrum is 



clearly dominated by emission lines. We also show a crude fit 
to the data using a non-equilibrium ionization model with vari- 
able abundances (model "vnei" in XSPEC). The fit is poor 
(X^ /v = 455 /254 = 1 .91), but this is mainly due to large resid- 
uals in the Ne emission line at 0.9 keV. In particular, the contin- 
uum component of the spectrum of RCW 89 is well-accounted 
for, with little emission above 3^ keV. 

While we defer a full and detailed treatment of these data 
to a subsequent paper, we think it clear from Figures ^ and 
Hthat we can rule out the possibility that RCW 89 consists 
of synchrotron-emitting clumps embedded in a diffuse thermal 
nebulae (as argued by G99), and that a model involving thermal 
clu mps embedded in a diffuse synchrotron nebula (as proposed 
by Tamura et al. 1996) seems far more likely. 



4. DISCUSSION 

4. 1 . The Overall Nebula 

It is immediately clear from the images presented here that 
the PWN surrounding PSR B 1509-58 has a very complicated 
morphology, matched only by the PWNe powered by the Crab 
and Vela pulsars in the level of structure. 

We first note that there is a clear axis of symmetry associ- 
ated with the system. This is manifested not only in the overall 
elongation of the nebula, but also in the orientation of features 
C, D and E in Figure ^, and features 1, 2, 4 and 5 in Figure ^ 
This axial structure is manifested on an extremely wide range 
of scales, from 10" = 0.2 pc up to 10' = 15 pc (where here and in 
further discussion we adopt a distance to the system of 5 kpc). 
It is hard to see how this orientation could be due to interaction 



with a pre-existing structure in the ISM (e.g. Caswell 1979) 



which should generally only influence the overall morphology 
of the system but not its small-scale structure. We therefore 
believe that this overall symmetry reflects the properties of the 
central pulsar. 

The only well-defined axes associated with a neutron star are 
its spin axis and its velocity vector. PWNe whose morphologies 
are dominated by the pulsar's motion are axially symmetric but 
inevitably s how a clear come t ary morphology wi th the pulsar 
at one end (^rail et al. 1996|; plbert et al. 2001|; K aspi et al. 
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2001). In contrast, the PWN surrounding PSR B 1509-58 ex- 
tends to many parsecs on both sides of the pulsar. It is thus 
clear that any motion of the pulsar cannot expla in th e overall 
PWN morphology (see additional discussion in §4.3), and that 
the main axis of the nebula must correspond to the pulsar spin 
axis. Such a correspondence was previously proposed for this 
pulsar by Brazier & Becker (1997), has al so been argued for 



the Crab and Vela p u lsars and their nebulae (Hester et al. 1995 
Pavlov et al. 2000| ; [Helfand, Gotthelf, & Halpern 2001D , and 



is seen in som e models for pulsar magnetospheres ( Sulkanen & 
Lovelace 1990). 

The region of the PWN where the symmetry of the system is 
most clearly broken is in the vicinity of feature F, the promon- 
tory extending several arcmin west of the main axis. The sharp 
edges seen for this feature are reminiscent of the "bay" ob- 
served alon g the western edge of the Crab Nebula (W eisskopf 
et al. 2000a). Such morphologies are presumably produced by 
strong confining pressure in these regions, possibly the result of 
pre-existing magnetic fields or circumstellar material. 

The spectral fits in Table || show that X-ray emission from 
the diffuse component of the PWN can be well fit by a power- 
law with photon index T = 2.05 ± 0.04 and absorbing col- 
umn (9.5 ± 0.3) X 10^' cm"^. This measurement is in good 



(iKawai, Okayasu, & Sekimoto 1993|; 


Irussoni et al. 1996 


Tamuraetal. 1996; Marsdenetal. 1997 


; Mineo et al. 2001). 



Du Plessis et al. (1995) have argued from a collection of 
archival X-ray data sets that the photon index for this PWN 
steepens from F « 1 .6 below e ~ 6 keV to F « 2. 15 above this 
energy. We see no evidence for this spectral curvature in our 
data set: a fit to the data using a broken power-law model rules 
out any spectral break across the Chandra bandpass larger than 
AF ^ 0.2. It is more likely that the effect claimed by du Plessis 
et al. (1995) was a result of the large uncertainties in and cali- 
bration differences between their lower sensitivity data sets. 

The standard observational picture for a PWN is that its spec- 
trum is comparatively flat at radio wavelengths, 1 < F < 1 .3, but 
is ste eper at X-ray energies, typically with F g:^ 2 (see G aensler 
2001 and references therein). It is usually assumed that syn- 
chrotron losses are at least partly responsible for this steepen- 
ing, generating a spectral break at energy Ej, across which we 
expect a change in ph oton index AF = 0.5 ( [Kardashev 1962 



Pacini & Salvati 1973). The frequency of this break, along with 
an estimate of the age of t he PWN, can be u sed to infer the 
nebular magnetic field (e.g. [Frail et al. 1996 ). In the case of 
the PWN around PSR B 1509-58, the observed photon index 
F = 2.05 implies Ej <; 1 keV. If the system is f yr old, then the 
implied nebular magnetic field is 



B„ = 8.0 



1700 yr 



-2/3 



V 1 keV / 



-1/3 



(1) 



so that we can infer 5„ ^ 8 /iG. 

There are a number of other, independent, estimates of the 
nebular magneti c field: the assumption of equipartition ( Se- 
ward et al. 1984), a simple model of PWN evolution (S eward 
et al. 1984) and a detection of inverse-Compt on emission from 
the PWN at TeV energies ( |Sako et al. 200C)| ) all result in esti- 
mates in the range B„ = 5-8 /iG, while Chevalier (2000) points 
out that the comparatively low efficiency with which this pulsar 
converts its spin-down energy into X-ray synchrotron emission 

"in fact the magnetic field in the inner parts of the PWN will generally be lower than the nebular average (§ [t.3.3| below; [Kennel & Coroniti 1984a| ), increasing the 
synchrotron lifetimes in these inner regions even further. 



also implies Ej 4 keV and hence B„ ;> 5 /iG. 

We thus conclude from a variety of arguments that the syn- 
chrotron break in this PWN is just below the X-ray band, and 
that the magnetic field is correspondingly low, Z?„ 8 /iG (cf. 
B„ ;> 100 /iG for the Crab). This low magnetic field pre- 
sumably results fro m the low confining pressure into which 
this PWN expands ( Seward et al. 1984 ). It is worth noting 
that if we assume a pulsar age f ^ 20 kyr (B landford & Ro- 
mani 1988; pvaramadze 200 1| ), we infer from Equation (|l]) that 
B„ ^1.5 /iG, which is at odds with all the other estimates de- 
scribed above. 

The above arguments imply that the synchrotron lifetime of 
electrons emitting in the Chandra band is similar to the age of 
the pulsar and nebula. Thus while synchrotron-emitting elec- 
trons in the outer diffuse nebula have now radiated most of their 
energy, we expect that regions closer to the pulsar will have a 
flow time less than the synchrotron lifetime and thus should 
have a spectrum flatter by AF ~ 0.5.''' Indeed the various X- 
ray features seen close to the pulsar (to be discussed in more 
detail in subsequent sections) all have photon indices F ^ 1.5. 

Figure |3] demonstrates a good match between the perimeter 
of the diffuse X-ray PWN and the "tongue" of emission appar- 
ent in radio observations. This correspondence was not appar- 
ent in earlier observations, which lacked the sensitivity to delin- 
eate the full extent of the X-ray PWN in this region (cf. Fig. 4 
of G99). Given that the two regions of emission occupy similar 
extents and have similar shapes, we propose that the "tongue" 
simply corresponds to emission from the pulsar nebula, and is 
the long-sought radio PWN in this system. This region has pre- 
viously been shown to have a very high fractional polarization 
and a well-ordered magnetic field ( Milne, Caswell, & Haynes 
1993; G99), properties typical of a radio PWN. 

The radio PWN is comparatively faint, sits on a complicated 
background, and is in close proximity to the very bright radio 
emission from RCW 89. The parameters of the radio PWN 
are thus poorly constrained, but using the data of Whiteoak 
& Green (1996) and of G99, we estimate an approximate flux 
density of 2 ± 1 Jy at both 0.8 GHz and 1.4 GHz. If we as- 
sume that this represents ^-^50% of the total radio flux den- 
sity of the PWN (the rest being hidden by RCW 89), we find 
/i GHz ^ (1 ±0.5) X 10"^ erg s"' cm"- keV"'. By comparing this 
to the flux density for the diffuse X-ray nebula inferred from Ta- 
ble H, fi keV = 1 .9 X 10"'' erg cm"^ keV"', we can infer that an 
unbroken power law with photon index F w 1 .6 is required to 
connect the two fluxes. This is is comparable to the photon in- 
dex for the inner emitting regions in the nebula which have not 
yet undergone synchrotron cooling, and is consistent with the 
spectrum 1 < F < 2 inferred for the radio emission ("region 6" 
in the analysis of G99). This good agreement with these other 
estimates of the uncooled synchrotron spectrum provides fur- 
ther evidence that the spectrum of the PWN must steepen from 
F ~ 1.5 to F « 2.05 just below the X-ray band, supporting the 
value B„ ~ 8 /xG inferred above. 

The radio observations shown in Figure |3] were designed to 
study the overall structure of the SNR, and were thus not op- 
timized to map out the fainter radio emission surrounding the 
pulsar Further radio observations are planned which we hope 
will allow us to better delineate the emission from this region. 

4.2. Outflow and Pulsar Orientation 
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Close to the pulsar, the collimated jet-like feature C is nar- 
row, bright, and closely aligned with the symmetry axis; at 
larger distances it fades, broadens and gradually changes direc- 
tion. Furthermore, Figure Q demonstrates that there is a deficit 
of radio emission along the length of the X-ray feature. These 
properties all argue that feature C corresponds a physical struc- 
ture rather than just an enhancement in the brightness of the 
nebula. We propose that feature C corresponds to a collimated 
outflow directed along the pulsar spin axis, as is thought to be 
seen for the C rab and Vela pulsars and PSR B0540-69 (H es- 
ter et al. 199 5; potthelf & Wang 20Qq; H elfand, Gotthelf, & 
Halpern 2001; p<^aaret et al. 2001 ), and as is also suggested by 
the morphology of the PWN in SNR G0.9-H0.1 ( Gaensler, Pivo- 
varoff, & Garmire 2001). The curvature of this outflow at large 
distances from the pulsar is reflected in the overall curved mor- 
phology of the surrounding nebula in this region (best seen in 
the ROSAT data of Fig. |l]). This curvature, as also seen in the jet 
in the Crab Nebula, presumably results from interaction with a 
gradient in the ambient density and/or magnetic field. 

In the case of the Vela pulsar, Radhakxishnan & Deshpande 
(2001) argue there is no true outflow along the spin axis, and 
that the apparent "jet" feature is a beaming effect produced by 
particles flowing along the pulsar's magnetic axis. However, 
in the case of PSR B1509-58 the observed curvature of the jet 
and the anti-correspondence with radio emission both provide 
good evidence for a physical cylindrical structure which, to per- 
sist over large spatial and temporal scales, can only be directed 
along the spin axis. 

As listed in Table ^, feature C has a photon index F « 1.7, 
distinctly flatter than the value F = 2.05 ± 0.04 determined for 
the diffuse surrounding component of the PWN. This contrast 
in spectrum indicates that the conditions in feature C must be 
distinct from those in the overall PWN: either the site and mech- 
anism for particle acceleration is different in the outflow from 
that in the diffuse PWN, or the synchrotron break in the outflow 
is at much higher energies. 

We first consider the possibility that there is continuous par- 
ticle acceleration along the outflow, which generates its harder 
spectrum. Such acceleration could be produced by shocks and 
turbulence in the interface between the outflow and its sur - 
roundings (Hjellming & Johnston 1988; [Frail et al. 1997 ). 
However, if the jet is optically thin (which we presume to be 
the case from its observed photon index) then we expect it to be 
brightest along its edges, not in its center as observed. We thus 
think this possibility unlikely. 

We therefore favor the alternative, that the main source of in- 
jection and acceleration in the outflow is the pulsar itself. There 
is no good understanding of how jets are produced in PWNe. 
For the purposes of this discussion, we assume that the injec- 
tion spectra for the jet and for the overall nebula are compara- 
ble. As discussed in §4.1, multiple lines of evidence indicate 
that the overall PWN cools from F 1.6 to F « 2.15 due to 
a synchrotron break at e,, <^ 1 keV. The spectrum for the jet is 
comparable to the uncooled spectrum for the diffuse PWN, in- 
dicating that any synchrotron break in the jet is at an energy 
£j > 10 keV ^ £j. We expect the magnetic field in the jet, Bj, 
to be significantly higher than that in the surrounding PWN, 
B„ ~ 8 fiG, in order to support the jet against surrounding pres- 
sure. Indeed we can estimate the minimum magnetic field in 
the outflow by assuming equipartition. We assume the jet to be 
a cylinder of length Ij = 4' = 5.8 pc and radius = 5" = 0.12 pc. 



corresponding to an emitting volume ~ 8 x 10 cm . From the 
spectral fit in Table |[ we infer an unabsorbed spectral lumi- 
nosity at 5 keV L5 keV = 1-3 x lO-'^ erg s"' keV"'. For a filHng 
factor (j> a nd a ratio of ion to electron energy /i, equipartition ar- 
guments ( Pacholczyk 197C| , ppl70-171) then imply a minimum 
energy in the jet Ej ^ 4-(l+ii)'*^''(jr'^'' x 10"^^ erg and a minimum 
magnetic field Bj - 25{l+ ^if/'^cj}-^/^ fiG > B„. The lifetime of 
synchrotron electrons at 5 keV is then tj < 14-0(1 + (jr"/^ yr, 
and so the rate at which the pulsar supplies energy to the jet is 
Ej ~ Ej/tj > 9(1 + fi) X 10^"* erg s"', coiTesponding to > 0.5% 
of the pulsar's total spin-down luminosity. 

The fact that the outflow shows an uncooled spectrum in 
the Chandra band implies ej ^ 10 keV. This can only be true 

if f/ < 100(1 + yr, where fy is the flow time in the 

jet. We can thus infer that the bulk velocity for the outflow is 
Vj = Ij/tf ^ 0.2c. The mean expansion velocity for the overall 
PWN (assuming free expansion into the SNR) is < 0.05c ^ vj. 
This confirms that the process by which the pulsar generates the 
jet is energetically and dynamically distinct from the process by 
which the overall pulsar wind is generated. 

The high jet speed we have inferred here is in stark contrast 
to the situation in t he Crab Nebula, whose polar outflow has 
a velocity < 0.03c dWilUngale et al. 200T1 ). While the nature 
of such jets is uncertain, we note that in the model of Sulka- 
nen & Lovelace (1990), a strong jet forms in the pulsar mag- 
netosphere only when the neutron star's magnetic field is suf- 
ficiently strong. We therefore hypothesize that the higher in- 
ferred dipole magnetic field for PSR B 1509-58 (a factor four 
higher than for the Crab pulsar) may be an explanation for its 
much more energetic outflow. The deficit of radio emission 
along the jet (seen in Fig. ^ is not easily explained; the particle 
spectrum in the jet may have a low-energy cutoff so t hat there 
are no radio-emitting pairs (e.g. Arons & Tavani 1993 ). 

The mean surface brightness"' for feature C in the energy 
range 0.3-8.0 keV is ~ 1 x 10"^ photons cm"^ sec"'. On the 
other side of the pulsar, in a region to the NW of feature E, the 
mean brightness is ~ 2 x 10"^ photons cm"^ sec"'. Any coun- 
terpart to feature C on the other side of the pulsar is therefore 
^ 5 times fainter. 

The absence of a counter-jet implies either that PSR B1509- 
58 is producing an intrinsically one-sided outflow, or that a 
counter-jet is present but cannot be seen. While multipole com- 
ponents to the n eutron star magnetic field could indeed produce 
a one-sidedjet ( Chagelishvili, Bodo, & Trussoni 1996 ), there is 
strong circumstantial evidence that the pulsar generates a colli- 
mated outflow to its north through which it is in teracting with 
the RCW 89 region ([Manchester & Durdin 19831; T amura et al. 
1996; prazier & Becker 1997| ; G99). This conclusion is fur- 
ther bolstered by feature D in Figure H, which is significantly 
elongated along the main axis of the system and widens with in- 
creasing distance from the pulsar. We propose that this feature 
represents a jet which produces no directly detectable emission, 
but for which we observe enhanced emission in a cylindrical 
sheath along the interface between the jet and its surroundings. 
Assuming that feature C and its unseen counterpart have similar 
intrinsic surface brightnesses and outflow velocities, relativistic 
Doppler boosting can account for the observed brightness con- 
trast provided that v^ cosf^ 0.28c, where ( is the inclination 
of the outflow to the line-of-sight. Since cos C < 1 for any in- 
clination, we can infer a lower limit Vj > 0.28c, consistent with 
the estimate vj > 0.2c determined above from the energetics of 



'^These values are absorbed surface brightnesses, but can be directly compared provided the absorption towards both regions is comparable. 



GAENSLER ET AL 



7 



the system. 

Brazier & Becker (1997) estimated > 70° on the basis of 
ROSATHRl data, in which the PWN appears to have a "cross"- 
shaped morphology suggestive of an edge-on torus. However, 
this inclination then results in an uncomfortably high velocity, 
Vj > 0.8c. In our Chandra observation, this cross-shaped mor- 
phology is not apparent; it presumably resulted from imaging 
of the inner and outer arcs (features E and 5 respectively) at the 
lower spatial resolution and sensitivity of the ROSAT HRI. A 



more reasonable post-shock flow speed 



:/3 (e.g. K ennel 

& Coroniti 1984a) corresponds to a smaller inclination angle 
( ~ 30°. This low value of ( has additional support from radio 
polarization observations of PSR B1509-58, which exclude the 
larger value s of C argued by Brazier & Becker ( 1997) at the 
~ 3(7 level ( [Crawford, Manchester, & Kaspi 200l[ ). 

Models in which the pulsar's high-energy emission origi- 
nates in outer gaps of the magnetosphere res ult in a viewing 
angle C ^ 45° (Romani & Yadigaroglu 1995), while those in 
which the 7-ray emission i s produced in a pola r cap generally 
require small values of ( (ECuiper et al. 1999). We are thus 
unable to distinguish between these models from the available 
data. 

4.3. Interpretation of Features E and 5 

Feature E, the bright semi-circular arc to the north of the pul- 
sar, demarcates a clear transition between the collection of com- 
pact bright features within 1' of the pulsar and the diffuse neb- 
ula which extends to much larger scales. The arc-like morphol- 
ogy of this feature is suggestive of a bow-shock, as would result 
where the ram-pressure from a fast-moving pulsar balances the 
outflow from the relativistic pulsar wind. For a bow-shock to 
result, the pulsar's motion must be supersonic with respect to 
the sound speed in the surrounding medium. Since the PWN 
itself has a sound speed c/ V3, this condition can only be met if 
the reverse shock from the surrounding SNR has collided with 
the P WN, bringing thermal material to the center o f the sys- 
tem ( IChevaHer 1998 ; [ van der Swaluw et al. 2()0lt B londin, 
Chevalier, & Frierson 2001). However, it is unlikely that this 
stage in PWN evolution has yet occurred in this system, since 
the resulting PWN would be brighter and more compact than is 
observed, and would occupy only a small fraction of the SNR's 
interior volume. 

An alternative interpretation is suggested by the fact that 
the dominant features seen in X-ray emission from both the 
Crab an d Vela PWNe are bri ght toroidal arcs surrounding the 
puls ar dHester et al. 1995|; H elfand, Gotthelf, & Halpern 
2001). In these cases, it is thought that these tori correspond to 
synchrotron-emitting particles from a pulsar wind focused into 
the equatorial plane of the system. Furthermore, in the case of 
the Crab, Chand ra data show a ring of emission situated inte- 
rior to the torus (Weisskopf et al. 2000a), which may represent 
the point where the free-flowing pulsar wind first shocks. 

In the case of PSR B 1509-58, features E and 5 both show an 
arc-like morphology which is bisected by the symmetry axis of 
the nebula, similar to what is seen in these other PWNe. The 
orientation inferred from the outflow being along the spin axis 
implies that features E and 5 are closer to the observer than is 
the pulsar, and are produced by a wind whose line-of-sight ve- 
locity component is directed towards us. The upper limit on 
any departure from circularity in the projected appearance of 
feature E implies C, <; 30°, in agreement with the estimate of 
the inclination angle determined in §4.2. Correcting for this 
inclination angle, we can infer a separation r^ = 0.4-0.5 pc be- 



tween the pulsar and feature 5, and = 0.75-0.85 pc between 
the pulsar and feature E. 

4.3.1. Analogy to the Crab Torus? 

In interpreting this emission, we first note that there are two 
characteristic time-scales associated with such features: tfiow, 
the time taken for particles to flow from the pulsar to this posi- 
tion, and tsynch, the synchrotron lifetime for these particles. 

We parametrize the upstream relativistic wind by the ra- 
tio, of the e nergy in electromag netic fields to that car- 
ried in particles ( Rees & Gunn 1974 ). This implies that a\ = 
B^/Airpijic^, where Bi is the toroidal magnetic field, pi is 
the lab frame total rest mass density, and cji ^ c is the four- 
velocity of the wind, all just upstream of the termina tion sh ock 
in the pairs. We assume cti ^ 1, and will show in § 4.3.3 that 
this assumption is self-consistent. In this case, the bulk velocity 
in the post-shock flow at a distance r from the pulsar is (K ennel 
& Coroniti 1984a): 



(2) 



where is the distance from the pulsar to the termination 
shock. We then find that: 



(3) 



For pairs emitting at energy e keV in a magnetic field B pG, 
the synchrotron lifetime is 



(4) 



For the Crab Nebula, we adopt rj sa 0. 15 pc and B 100 pG, 
and consider the Crab's X-ray torus at r w 0.4 pc. At an energy 
e = 5 keV, we find f//o„, ~ 30 yr and tsy„ch ~ 20 yr, so that these 
time scales are comparable. In the case of PSR B 1509-58, we 
assume that the location of the termination shock corresponds 
approximately to that of feature 5 (the inner arc), and so adopt 

« rg . For a magnetic field B ~ B„ = 8 fiG, we then find that in 
feature E, tfiow ~ 25 yr <C fsv«c/i ^ 800 yr Since the post-shock 
magnetic field in the inne r nebula is generally weak er than the 
mean value for the PWN (Kennel & Coroniti 1984a), this value 



of fjync/i is likely to be a lower limit, further widening the dis- 
crepancy between the two time-scales. 

Thus while the bright torus seen in the Crab simply corre- 
sponds to the region in which most of the X-ray emitting parti- 
cles radiate their energy, synchrotron cooling is not significant 
at this distance from PSR B1509-58 even at X-ray energies, 
and the brightness enhancement in feature E cannot be due to 
rapid dumping of pairs' energy into X-ray photons. The long 
synchrotron lifetimes require that most of the X-ray emission 
comes from a much larger volume, as is indeed observed. 

It is therefore not valid to argue that feature E is the analog 
of the X-ray torus seen in the Crab; the much lower nebular 
magnetic field in the case of PSR B1509-58 demands a differ- 
ent interpretation. This conclusion is supported by the fact that 
feature E has a distinctly harder photon index than the overall 
nebula and also shows a clear radio counterpart (see Fig. Q), nei- 
ther of which would be expected if this feature were primarily 
due to radiative losses. 
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4.3.2. "Wisps" in an Equatorial Flow 

The Crab Nebula shows a series of toroidal optical and X-ray 
features cl ose to the pulsar, generally termed "wisps" (L am- 
plan d 1921~|Scargle 1969|; [Hester et al. 1995|; M ori et al. 
2002). These wisps form at the termination shock, move out- 
wards in the equatorial plane at ~ 0.5c, and then slow and fade 



netosphere's potential is 



(Fanvir, Thomson, & Tsikarishvili 1997; Hester et al. 1996) 



Furthermore, radio wisps have recently been identified in the 
Crab Nebula, originating in the same region as the wisps seen 
at high er energies, and similarlv moving outwards at high ve- 
locity ( Bietenholz, Frail, & Hester 2001 ). This result suggests 
that the radio, optical and X-ray emitting populations are all 
accelerated in the same region, a result difficult to explain in 
all models advanced to date for shock acceleration at a PWN's 
wind termination shock. 

As we have shown in § [4.3. 1| , the much lower magnetic field 
for the PWN around PSR B1509-58 results in hynch > f//on for 
feature E, implying that it cannot be a clone of the Crab's torus. 
In the following discussion, we therefore explore the conse- 
quences of identifying features 5 and E as analogs of the wisps 
seen in the equatorial outflow from the Crab. The fact that we 
see in Figure a clear radio counterpart to X-ray feature E, as is 
seen for the wisps in the Crab, gives phenomenological support 
to its interpretation as a wisp-like feature, while the high level 
of linear polarization is consistent with the expected toroidal 
geometry for the magnetic field in this region. We note that be- 
cause the synchrotron cooling time is long compared to the flow 
time, thermal instabilities due to synchrotron cooling, proposed 
by Hester (1998) to account for the wisps in the Crab Nebula, 
cannot account for features E and 5 seen here. 

In the Crab Nebula, Gallant & Arons (1994, hereafter 
GA94) extended the work of Kennel & Coroniti (1984a) by 
modeling the first optical wisp as emission from magnetized 
electron/positron pairs, compressed by magnetically reflected 
heavy ions just behind the leading edge termination shock in the 
pairs. These ions form a highly accelerated Goldreich- Julian 
return current in the flow. For such magnetic reflection in the 
overall shock structure to be a viable explanation of the wisps 
observed in the Crab, the ions must have energy per particle 
comparable to the maximum energy possible, 7,m,c^ w Ze^open 
(where 7,-, m, and Ze are the Lorentz factor, mass and charge of 
the ions respectively, and $open is the open field potential of 
the pulsar's magnetosphere). S ubsequent time-dependent stud- 
ies ( Spitkovsky & Arons 2000 ) have shown the time-variability 
due to the electromagnetic instability of the reflected ion cur- 
rent to be remarkably similar to t he temporal behav ior reported 
for the Crab's inner X-ray ring ( Mori et al. 2002 ) and possi- 
bly f or the first optical wisp (T anvir, Thomson, & Tsikarishvili 
1997; fester et al. 1996| ). 

In the following sections, we propose that feature E is an ana- 
log of the Crab's second wisp, and results from compression at 
the second, disorganized ion turning point in the pulsar's equa- 
torial outflow. 

4.3.3. Structure of the Flow 
In the model of GA94, the upstream flow Lorentz factor is: 



71 = ri — = 2.6 X 10 



0.33 A Vl-8x 10" erg s- 



1/2 



(5) 

where 77 is the fraction of the open field potential picked up by 
the ions, A is the atomic number of the ions, and the open mag- 



^open ~ 



1/2 



:7.4x 10'" V. 



(6) 



Previous uses of this model to study X-ray emission from the 
Crab, PSR B1957H-20 and PSR B 125 9-63 have consistently 
found ?7 ^ 1/3 to be the bes t choice (A.i'ons & Tavani 1993; 



GA94; [Tavani & Arons 1997[ ). Soft X-ray observations suggest 
that hydrogen or possibly helium domi nate the composition o f 
the upper layers of pulsar atmospheres ( Pavlov & Zavlin 2000 ), 
from which the ions of the return current are drawn. We conse- 
quently adopt Z ^ 1 and A 1 in further discussion. 

The intrinsic thickness of the pair shock is comparable to the 
Larmor radius in the pairs, TZ\±, based on the upstream Lorentz 
factor, 7i, and upstream magnetic field, B\ = B2/3, where B2 is 
t he magnetic field im mediately downstream of the pair shock 
( [Gallant et al. 1992[ ). Adopting B2 = 2.5 ^G (see Equation [|] 
below) and 71 - 2.6 x 10*, we find that 7^l± 5 x 10"^ pc, 
corresponding to an angular scale at 5 kpc of 0."03. The pair 
shock crea tes a relativistic thermal distribution of pairs (G allant 
et al. 1992), with mean energy per particle ^ ^inieC^ ^ 1 TeV, 
whose synchrotron emission in the immediate post-shock mag- 
netic field is in the infrared, peaking around 5 /im and unob- 
servably faint because of the slow cooling. The observed non- 
thermal X-ray emission from feature 5 and at larger radii is due 
to the acceleration of pairs from the thermal pool. Therefore, 
in this interpretation, feature 5 is the analog of the first wisp in 
the Crab Nebula, and is created by the deposition of the heavy 
ions' outflow momentum at the first turning point in their mag- 
netically reflected orbits. 

Identifying feature 5 as an ion-driven compression implies 



1 1 m,c-7i 

'"^-2^^'=2 Z^ 



1 V^open 

2 B2 



where 7?,2, is the downstream ion Larmor radius. 
r^ ~ 0.5 pc, we then find that 



3Bi=B2 = 



open 



2r5 



= 2.5^ 



0.33 



/iG. 



(7) 



Adopting 



(8) 



With a magnetic field at feature 5 of 2.5 fiG, the Lorentz 
factor of the X-ray emitting pairs is 7 w 4 x 10^, correspond- 
ing to a Larmor radius TZ2± ~ 0.09 pc. Therefore, the finite 
gyration of the pairs emitting in X-rays will cause feature 5 to 
have an angular width ^ 4" at an energy of 5 keV, no mat- 
ter what the intrinsic thickness of the structure in the magnetic 
field and the pair plasma. Figure || shows that feature 5 is in- 
deed 4" -5" across. Most of the pairs in the flow have energy 
well below the energies of the X-ray emitting particles. In a dis- 
tribution for which the number of particles N±{j) cx 7"'', where 
/:> = 2r - 1 ^ 2 these lower-energy particles are just more effi- 
ciently radiating tracers with no substantial dynamical signifi- 
cance. 

The radius of the next compression, r„ex/, can be readily es- 
timated. We expect B(r„exi) = kB{rs){rnext/r5), where the factor 
^ > 1 represents the ion induced compression of the magnetic 
field above the linear increase of B with r; the linear increase 
occurs simply because of deceleration of the expanding flow. 
GA94 found k ^ 2, while the time dependent simulations of 
Spitkovsky & Arons (2000) show /t ~ 1 .5. 
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The radius of the second compression occurs approximately 



at 



r5 + TZi[B(r„ext)] ='"5 + 



next 



(9) 



Adopting ^ ~ 1.5 yields r„e„/rs = 1.76, or r„„f = 
0.88 (rs/O.S pc) pc. The good agreement between r„„, and 
the observed radius for feature E, = 0.75 - 0.85 pc, sup- 
ports identifying feature E as the second ion compression in 
t he flow. Relativistic cyclotr on instability of the reflected ions 
( [Spitkovsky & Arons 2000| and in preparation) disorganizes 
their flow at larger radii, and no further coherent compressions 
(a.k.a "wisps") at still larger radii can be produced by this 
model, in agreement with the observations. 

The model predicts the peak magnetic field in feature E to be 



BE = kB2 — 



6.5 



0.5pc T] 
rs 033 



/zG; 



(10) 



the average magnetic field in feature E is closer to 5 /iG. 



from the star flow out along field lines which fill the whole 
space around the star, in the polar as well as the equatorial sec- 
tors. Therefore, the pairs expand quasi-spherically so that 

A57 K± _■, 

ni± = «!+ + «!_ = 2K±nu- — = 6.9 x 10 cm , (15) 

Att 10^ 



and so 



^(theoretical) ^ __A__ ^ 4 ^ IQ'^ — ( _±_ ) ; (16) 



0.33 / rs 



47rn i± 71 m^c^ 



77 \0.5pc 



the ions have five times the flow energy of the pairs. 
Since 

1 _ 1 ^ 1 

0"! ai± fJi,' 



we find that 



^(theoretical) _ y ^ -|^Q-4 V ''5 



0.33 0.5 pc' 



(17) 



4.3.4. Energetics of the Flow 

The energetics of the ion flow can be estimated from theory, 
since the electrodynamics of the magnetosphere constrain their 
outflow rate to the Goldreich- Julian value. 



Ni 



2^openC 

Ze 



: 3 X 10"" s 



33 „-l 



From Equation (yj), the ion density just upstream of the shock 
in pairs'^ is then 



where we have assumed k± = 10^. 

Estimating the pair density and thus ai± from X-ray data 
is more difficult, since the relatively narrow bandwidth of the 
observations means that we have very little constraint on the 
overall energetics. The spectral fits to feature E listed in Ta- 
ble H indicate an unabsorbed 0.5-10 keV flux density /e ~ 
(11) 3.6x10-'^ erg cm s and a photon index F = 1 .7. For a dis- 
tance to the source of 5 kpc, these imply a spectral luminosity 
in the Chandra band of 



nii 



= 3.3 X 10 -TTT cm 

An 



Airr^c An 



(12) 



where Aft is the solid angle subtended by the ion outflow. In 
the Crab, the appearance of the wisps and of the X-ray torus 
suggest 4tt/ An ^ 5.8. We do not know whether the equatorial 
outflow from PSR B1509-58 is similarly flattened, although the 
presence of the jet (feature C) clearly indicates that some degree 
of flattening is present. 

The most meaningful way to characterize the ions' energet- 
ics is to compare them with the electromagnetic energy flux just 
upstream of the pair shock, by computing the parameter a for 
the ions. We find 



(theoretical) _ 
ll 



47rni/7im/c^ 



: 8 X 10" 



0.33 5.8 



4tt_ 

An 



0.5 pc 



(13) 

Theory suggests that PSR B 1509-58 should produce k± w 
10^ electron/positron pairs per primary particle in the polar 
current-carrying beam within the magnetosphere if the surfac e 
magnetic field is strictly dipolar ( ^ibschman & Arons 2001 ); 
K± might be a fe w times larger if the sur face magnetic field is 
stron gly distorted (|Barnard & Arons 1982|; Hibschm an & Arons 
2001) or if the dipole is strongly offset ( [Arons 1998 ). Then 



y(theoretical) 



: K±Ni : 



:3xlO-^s- 



:2.8x 10^^ e-°-^ergs-' keV"'. 



(18) 



The power law form of the spectrum of feature E suggests 
synchrotron radiation from a power law distribution of parti- 
cles. From synchrotron theory, the isotropic emission rate of an 
isotropic power law distribution of electrons and positrons with 
ener gy spectrum N+(E) = KE~'\ E,j,i„ <E < E^ax is (H jellming 
1999) 



Is 



1.1 X IQ-^'^ KB],'^ e'^-^ ergcm"^ s"' keV"', 



(19) 



where we have used p = 2r - 1 = 2.4. 

We model feature E as occupying a volume Ve = 
47rr|Ar(AJ7/47r), where Ar is the radial width of feature E. 
Adopting r£ = 0.8 pc, Ar w 0.4 pc and 47r/AO = 5.8, we find 

= 1.6 X lO^*" cm-'. Since = VEqe, we can combine Equa- 
tions (|l|) and (|l|) to find that ^ = 1 .0 x 10"^ erg' ^ cm-^ where 
we have used Be = 5 /iG. Integrating the power law particle 
spectrum over energy yields a total pair density in feature E of 



+ «i^' = 2.3x lO^El;'*! 



(20) 



where E^m eV is the minimum energy of the power-law particle 
distribution. The injection rate into feature E is therefore 



N^E) ^ 4^^2^„(£)^(^^) ^ J 2 X 10"^ E^ s 



^54 



1.4 „-l 



(21) 



(14) As per Equation (||), we used v(r£) = (c/3)(rs/rE)^ = 3.9 x 



The ion compressions that lead to the observed wisps only 
occupy an equatorial sector of latitudinal extent AO where the 
electric return current from the pulsar flows, while the pairs 

"Assuming rj is not much less than the radius of the first turning point in the ions' 



10^ cm s"' in the numerical evaluation 

We have no observational handle on /Smin, other than it 
must be small enough to allow the synchrotron critical en- 
ergy (3/2)(heB/mec)(Emin/'neC^)^ to be well below the lower 

orbits, as is true for the GA94 model of the Crab. 
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limit of the effective Chandra band for this source, e ~ 1 keV. 
Therefore, observationally < 53 TeV. Theoretically, ex- 
isting work on the particle acceleration properties of quasi- 
transverse relativistic magnetosonic shocks suggests the down- 
stream power law shoul d flatten (and perhaps cu t off) be- 
low the en ergy ^im^c ^ (iKennel & Coroniti 1984b|; Hoshino 
et al. 1992). From Equation (|^), we then would have Ej^in > 
1.3(77/0.33) TeV. The particle acceleration theory suggests the 
true power form of the particle distribution sets in at two to 
three times this energy, in which case 



'4x10 



,36 



0.33 



(22) 



which compares well with Equation (|l4|). 
Equation (E3) yields: 



; 7 X 10" 



033 



(23) 



a similar value to the theoretical estimate of Equation ([Tq). 
From either approach (Ti± is small, and is comparable to the 
value (T w 3 x 10"^ inferred for the Crab Nebula (K ennel & 
Coroniti 1984a7^ mmering & Chevalier 19871). 

Feature E forms less than half a complete torus — this arc 
is 3-5 times brighter than any counterpart on the other side of 
the pulsar The flow velocity given by Equation (||) is too small 
at feature E to invo ke Doppler boosting as an explanation (cf. 



Felling et al. 1987). We thus can offer no explanation for this 



asymmetry, except to note that a similar situation is seen in the 
equatorial sti'uctures in the Vela PW N and in SNR G0.9H-0.1 
( [Helfand, Gotthelf, & Halpern 200 1|; G aensler, Pivovaroff, & 
Garmire2001). 

Nevertheless, we believe that the ideas developed by GA94 
and Spitkovsky & Arons (2000) to explain the spatial (and tem- 
poral) structure of the equatorial interaction zone between the 
Crab Nebula and its pulsar can be usefully applied to the equa- 
torial interaction between PSR B 1509-58 and its PWN, as re- 
vealed by the Chandra observations presented in this paper. To 
the extent that the model outlined here does represent the facts, 
the basic conclusion is the same as was the case in the GA94 
model for the Crab, namely cti <C 1 . This result has been ex- 
tracted from a close study of the detailed structure of the in- 
teraction zone rather th an from modeling the glo bal dynamics 
of the whole PWN (cf. fennel & Coroniti 1984a| ), the latter of 
which is difficult to apply to the complicated morphology sur- 
rounding PSR B 1509-58. 

We note that if features E and 5 are indeed analogs of the 
wisps seen in the Crab, then we similarly expect them to move 
away from the pulsar at high velocity. For example, motion 
outwards at 0.5c would correspond to proper motions of a few 
arcsec per year at the distance of PSR B1509-58, which could 
easily be detected by Chandra at subsequent epochs. 

4.4. Small Scale Structure Near the Pulsar 

The high resolution of Chandra has revealed emission from 
several small-scale features close to the pulsar, as seen in Fig- 
ure | 

Feature 6 is likely to cor re spond to emission from t he O star 
Muzzio 10 dMuzzio 19791; prsatti & Muzzio 1980|). X-ray 



by a Raymond-Smith spectrum with kT w 0.2-0.5 keV and an 
unabsorbed luminosity (0.3-10.0 keV) of ^ 1 - 3 x 10^^ erg s"' 
(Berghofer, Schmitt, & Cassinelli 1996), corresponding to 
an unabsorbed flux density (0.3-10.0 keV) fx ^ (0.4-4) x 
lO"'"* er g cm~^ s~' for a pho t ometric distan ce in the range 2.5- 



4.6 kpc ( jSewai-d et al. 1983| ; [Arendt 199 1| ). The crude spectral 
parameters inferred for this source in Table || are consistent with 
these values. 

What features 1^ represent is no t immediately clear. If we 
accept the argument made in §4.3.1 that rj w rs, then features 
1-3 (and possibly feature 4, depending on projection effects) 
originate in the zone in which the wind is still freely expand- 
ing. 

In the Crab Nebula, a variety of small-scale structures has 
similarly been identified within the unshocked wind zone. 
"Knot 1" and "knot 2" (the latter of which is also termed "the 
sprite") are resolved optical structures lying 1500 and 9000 AU, 
respectivel y, from the Crab pulsar along the jet axis ( Hester 
et al. 1995). It has been proposed that these features corre- 
spond t o quasi-stationary s hocks in t he polar outflow from the 
pulsar ( [Vlichel et al. 1996 ; Lou 1998), or, in the case of knot 1, 
to a sheath of emission surrounding this outflow (G raham et al. 
1996). The two knots both show significant variability in their 
brightness, position and morphology on tim e scales of days 
to months ( j^ester et al. 1996 ; Hester 1998). Similarly time- 
variable kno ts of emission are also seen in X-rays close to the 
Vela pulsar ( |Pavlov et al. 200 1| ). 

In the following discussion, we focus on feature 1, the knot 
of emission closest to PSR B1509-58. We estimate feature 1 
to have approximate dimensions 3" x 5", and its projected sep- 
aration from the pulsar to be r^not ~ 4" =0.1 pc. We assume 
that this emission is generated by particles accelerated in some 
localized turbulent region, such as might be produced by the 
collision of inhomogeneous wind streams proposed by Lou 
(1998). The spectrum observed for feature 1 (F ~ 1 .2) is some- 
what harder than the uncooled spectrum for the overall PWN 
(F ^ 1 .6), supporting the possibility that this is a region of sep- 
arate particle acceleration. In this case, we can interpret the 
extent of this feature as corresponding to the relativistic Lar- 
mor radius of gyrating pairs. Although the number of photons 
available is limited, there is the suggestion in the data that the 
extent of the knot increases slightly with increasing photon en- 
ergy, consistent with this interpretation. At a photon energy 
e keV and in a magnetic field Bknot the pair Larmor radius 
is 7^knot = 0.13(e/B3^ji/2 pc. Adopting 7^knot 0.05 pc and 
e = 5 keV, we can in fer fik not ^ 3 fiG. 

We computed in §4.3.4 that the rate of pair production in the 



pulsar wind is N± « 4 x 10 s . The pair density in feature 1 
is therefore 



W± knot — 



N± 



1.1 X 10" 



■'«cm-^ 



(24) 



emission from an 06.5111 star can be typically approximated 



In Equation (|]) we computed an upstream flow Lorentz fac- 
tor 71 = 2.6 X 10^. Combining these estimates for Bknot, n±.knot 
and 71, we can infer by analogy with Equation ( [l6| ) that a± < 
3 X lO""* in this region. (We adopt this value as an upper limit, 
since compression and turbulence likely enhance the magnetic 
field in this feature above the ambient value.) 

This low value of a is consistent with the various models for 
energy transport in the unshocked wind, which generally re- 
quire the transition from cr ^ 1 (at the light cyli nder) to cr <C 1 
(at the termination shock) to occur at r <C ^knot (Coroniti 1990 
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Michel 1994 ; Melatos, private communication). It is interest- 
ing to note that in the large-amplitude plasma wave model pro- 
posed by Melatos (1998), the wind beyond this transition point 
evolves as cr cx r^. We thus expect a for feature 1 to be ~20 
times smaller than that inferred at the termination shock, a re- 
sult not inconsistent with our calculations. 

A better insight into the nature of these compact features 
will require observations of this system at multiple epochs and 
at longer wavelengths (e.g. in the near-infrared). With such 
data, we will be better able to determine the broadband spec- 
trum and hence total energetics of these sources, can establish 
whether these features are persistent or transient, and can de- 
termine whether any of these sources shows outward proper 
motion which could associate them directly with outflow from 
the pulsar Indeed, Chandra has identified significant changes 
in the positions and brightnesses of small-s cale structure in th e 
Vela PWN on a time-scale of seven months (Pavlov et al. 2001 ). 



5. SUMMARY AND CONCLUSIONS 

Our Chandra observations have confirmed that PSR B1509- 
58 interacts with its surroundings in a spectacular and highly 
anisotropic fashion. Our main results are as follows: 

1. PSR B1509-58 powers a large elongated non-thermal 
nebula, whose main symmetry axis most Ukely corre- 
sponds to the pulsar spin axis. A variety of arguments 
implies a mean nebular magnetic field of B„ ^ 8 /iG, 
corresponding to a synchrotron break just below the X- 
ray band. The high sensitivity of these data has allowed 
us to map out the full extent of the X-ray PWN. We 
have shown that this nebula matches well the extent of 
a previously-identified "tongue" of radio emission near 
the pulsar, and propose this latter feature to be the radio 
nebula powered by pulsar 

2. We confirm previous claims of a one-sided coUimated 
outflow emanating from the pulsar along the main sym- 
metry axis, and which generates a channel of reduced 
radio emission. We determine that this jet carries away 
at least 0.5% of the pulsar's spin-down luminosity, and 
that the bulk velocity in the outflow is v j > 0.2c. We 
have identified an elongated sheath on the other side of 
the pulsar which indicates the existence of a counter-jet, 
despite the lack of direct emission from such a source. 
The brightness contrast between the two outflows can be 
explained by Doppler boosting, and implies an inclina- 
tion between the pulsar spin axis and the line-of-sight 
C--30°. 

3. Approximately 0.5 and 0.8 pc to the north of the pulsar, 
we identify two arcs of emission with a toroidal geome- 
try which similarly implies an inclination ( ^ 30°. A re- 
analysis of radio observations of the region shows that the 
second of these arcs has a clear radio counterpart and is 
strongly linearly polarized. We show that the flow time in 
these features is much less than the time scale for radia- 
tive losses, and that they therefore cannot correspond to 
a radiating torus as seen in the Crab Nebula. Rather, we 



interpret these features as "wisps" in an equatorial flow, 
and show that such an assumption is consistent with these 
sources representing internal structure in the termination 
shock of a particle-dominated ion-loaded wind, as pro- 
posed by Gallant & Arons (1994) to explain the wisps in 
the Crab Nebula. We infer a ratio at the shock of Poynting 
to particle flux for electron/positron pairs cr± ~ 0.005, a 
value similar to that argued for the Crab. 

4. We identify several compact knots of emission less than 
0.5 pc from the pulsar, within the unshocked wind zone. 
The hard spectrum seen for these knots suggests a differ- 
ent acceleration mechanism than for larger scale nebular 
features beyond the termination shock. We infer a ratio 
of Poynting to particle flux for one of these knots a± < 
0.003, indicating that the transition from a magnetically- 
dominated to a particle-dominated wind occurs less than 
0. 1 pc from the pulsar 

As has been the case in the past, new data on this source raise 
as many questions as they answer. We lack an understanding of 
how PSR B1509-58 generates such a fast-moving and colli- 
mated outflow over such large scales, nor is it clear why such a 
striking feature is not seen in most other PWNe. We have also 
not addressed the nature of the interaction between this outflow 
and the surrounding SNR; a forthcoming paper will present an 
analysis of X-ray emission from the RCW 89 region, which 
may better elucidate the physical conditions associated with 
this process. While we have outlined a self-consistent model 
in which the arcs seen near the pulsar are wisps as in the Crab, 
we have no explanation for why these features are only seen 
on one side of the pulsar (it is important to realize that the na- 
ture of the Crab's wisps are themselves still a matter of heated 
debate). We clearly need to extend our coverage to other wave- 
lengths and additional epochs if we are to better understand the 
observed structures. Finally, the process which produces the 
innermost X-ray features is completely unknown, and will also 
require further study. 

To conclude, we comment that for many decades the Crab 
Nebula has been the focus of almost all efforts to understand 
the process by which pulsars couple to their environment. With 
data such as those presented here, we can now finally extend 
such studies to a wider sample of sources. 
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Table 1 

Detected count rates for selected regions. Uncertainties are ±la. 



Region 


Count Rate"" 
(cts s"') 


Mean Surface Brightness'' 
(10"^'' cts s"' arcsec"^) 


Entire source 


8.79 ±0.03 


2.04 ±0.01 


PSRB 1509-58 


0.140 ±0.003'' 




RCW 89 


2.39±0.01 


3.51 ±0.01 


C 


0.212 ±0.005 


3.38 ±0.08 


E 


0.146 ±0.005 


3.6±0.1 


1 


0.014±0.001 


78±6 


2 


0.007 ±0.001 


29 ±4 


3 


0.042 ± 0.002 


86±4 


4 


0.029 ±0.002 


66±5 


5 


0.068 ±0.003 


15±1 


6 


0.0017 ±0.0004 





"Count rates and surface brightnesses are those detected over the energy range 0.5-10 keV. and have been corrected for background. Surface brightnesses are only 
given for extended sources. 

''The pulsar suffers from severe pile-up; this value is therefore a significant underestimate of the incident count rate. 



Table 2 



Spectral fits to various sub-regions of the source. 



Region 


Model 


A^H (10-' cm"^) 


r / kT (keV) 


/, (10"'^ ergcm"^ s"')'' 




Diffuse PWN 


PL 


9.5±0.3 


2.05 ±0.04 


55 ±3 


1686/1268= 1.33 


C 


PL 


8±1 


1.6±0.1 


5±1 


212/237 = 0.89 




PL 


9.5 (fixed) 


1.70 ±0.09 


5.2±0.5 


217/238 = 0.91 


E 


PL 


8!f 


1.6 ±0.2 


3±1 


206/193= 1.07 




PL 


9.5 (fixed) 


1.7±0.1 


3.6 ±0.4 


207/194= 1.07 


1 


PL 
PL 


10!^ 
9.5 (fixed) 


^•■'-0.5 

1.2±0.3 


'-'•^-0.4 

0.4±0.1 


24/21 = 1.13 
24/22= 1.08 


2 


PL 


9+20 


6+1 
^•"-0.6 


-0.2 


12/11 = 1.10 




PL 


9.5 (fixed) 




-0.2 


13/12= 1.06 


3 


PL 
PL 


11+3 

9.5 (fixed) 


1.7±0.3 
1.6±0.1 


^ -^-0.3 

1.0 ±0.2 


57/61 = 0.93 
57/62 = 0.92 


4 


PL 
PL 


iii 

9.5 (fixed) 


1.3 ±0.4 
1.3±0.2 


0.8±0.2 


50/52 = 0.96 
50/53 = 0.95 


5 


PL 
PL 


6±2 
9.5 (fixed) 


^■■'-0.2 

1.6±0.1 


^•"-0.5 

1.6 ±0.2 


101/84= 1.21 
105/85= 1.24 


6 


RS 


-10 


-0.6 


-0.1 





Uncertainties are all at 90% confidence. All models assume interstellar absorption using the cross-sections of Baluciiiska-Church & McCammon (1992), assuming 
solar ahiindanr ps Mndpis used- "PI " i nriiratps a power law of the form /e oc e"^ where r is the photon index; "RS" indicates a Raymond-Smith spectrum of 
temperature T (llaymond & Smith 19771. 

"Flux densities are for the energy range 0.5-10 keV, and have been corrected for interstellar absorption. 
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Fig. 2. — Chandra image of G320.4-1.2 over the entire energy band. The image has been exposure-corrected and convolved with a gaussian of FWHM 10" 
and is displayed using a linear transfer function. The central black box in the main image shows the region covered by the panel at lower right. In this panel, thi 
data have been smoothed with a gaussian of FWHM 4" and displayed using a logarithmic scale; the white box in this panel shows the region displayed in Figure p 
Specific features discussed in the text are indicated. 
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Fig. 3. — A compai'ison of X-ray and radio emission in tlie region surrounding PSR B1509-58. Tlie greyscale displays radio data from 1.4-GHz observations 
with the ATCA at a spatial resolution of 21" X 24" and over a range to 10 mJy beam"' , as presented by G99. The tipples seen in the radio emission are residual 
sidelobes from the bright RCW 89 region to the north. The contours show Chandra data as in Figure O, with contours at 0.5%, 2% and 30% of the peak. 
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Fig. 5. — X-ray emission immediately surrounding PSR B 1509-58 in the energy range 0.3-8.0 keV. These data have been exposure-corrected but have had no 
smoothing applied. Specific features discussed in the text are marked and labeled. 
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Fig. 7. — A low-resolution version of Figure ^ showing the extraction regions used to generate the spectra shown in Figures ^ and ^. The region marked "Diffuse 
PWN" is an annular region which excludes regions A, C and E and features 1-5. 
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Fig. 8. — X-ray spectrum of the diffuse pulsar wind nebula. The data-points in the upper panel show the data from each of the four CCDs on which the nebular 
emission falls, while the solid lines correspond to the best-fit absorbed power-law model. The lower panel shows the number of standard deviations by which the 
model and the data differ in each bin. The data have been plotted so as to give a signal-to-noise ratio of 10 in each bin. 
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Fig. 9. — As in Figure ^ but for a spectral fit to the brightest clumps in RCW 89. In this case, the sohd line corresponds to a non-equilibrium ionization model 
with variable abundances, and the data have been plotted with a signal-to-noise of 8 in each bin. 



